Oblique surface Josephson plasma waves in layered superconductors 
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We have theoretically studied oblique surface waves (OSWs) which propagate along the interface 
between a dielectric and a layered superconductor. We assume that this interface is perpendicular 
to the superconducting layers, and OSWs at the interface can propagate at an arbitrary angle 
with respect to them. The electromagnetic field of the OSWs in a layered superconductor is a 
superposition of an ordinary wave (with its electric field parallel to the layers) and an extraordinary 
wave (with its magnetic field parallel to the layers). We have derived the dispersion equation for 
the OSWs and shown that the dispersion curves have end-points where the extraordinary mode 
transforms from evanescent wave to bulk wave, propagating deep into the superconductor. In 
addition, we have analytically solved the problem of the resonance excitation of the OSWs by the 
attenuated-total-reflection method using an additional dielectric prism. Due to the strong current 
anisotropy in the boundary of the superconductor, the excitation of the OSWs is accompanied by 
an additional important phenomenon: the electromagnetic field component with the orthogonal 
polarization appears in the wave reflected from the bottom of the prism. We show that, for definite 
optimal combinations of the problem parameters (the wave frequency, the direction of the incident 
wave vector, the thickness of the gap between dielectric prism and superconductor, etc.), there is 
a complete suppression of the reflected wave with its polarization coinciding with the polarization 
of the incident wave. Contrary to the isotropic case, this phenomenon can be observed even in the 
dissipationless limit. In such a regime, the complete transformation of the incident wave into a 
reflected wave with orthogonal polarization can be observed. 

PACS numbers: 74.78. Fk, 74.50.+r, 74.72.-h 



I. INTRODUCTION 

Surface electromagnetic waves represent a specific kind 
of macroscopic perturbations which propagate along in- 
terfaces between different media. The electromagnetic 
field in such waves decay exponentially away from the 
interface deep into both media. A very important exam- 
ple of surface waves is the so-called plasmon-polaritons 
observed at the interfaces between a normal metal and a 
dielectric at frequencies about the plasma frequency uj p , 
which belongs to the optical or far-infrared ranges 1 - - — . 
Plasmon-polaritons are central in numerous resonance 
phenomena. Examples of these are: (i) the extraordinary 
transmission of light through metal films with subwave- 
lcngth holes^; (ii) an "inverse" effect of the resonant 
suppression of light transmission through modulated 
metal films with thicknesses less than the skin-depth 7 ; 
(iii) the Wood anomalies in the reflectivity^— £ and trans- 
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of periodically-corrugated metal sam- 



ples. These optical anomalies could have potential appli- 
cations for photovoltaics, light control, filtering and de- 
tection of radiation in far-infrared and visible frequency 
ranges. 

It would be very desirable to similarly control the elec- 
tromagnetic radiation in the terahertz (THz) frequency 
range. The mastering of this range (0.3-10 THz) is a 
rapidly developing area of research due to promising ap- 
plications. In principle, plasmon-polaritons can exist at 
w « w p . However, the dispersion curves of plasmon- 



polaritons are very close to the light-line uj = ck/ -Je2 in 
the THz frequency range. In this case, the surface waves 
are almost extended^, in the sense that they are very 
weakly localized in the dielectric over distances of about 
1 meter (here u and k are the frequency and wave vector 
of plasmon-polaritons, Ed is the permittivity of the dielec- 
tric, and c is the speed of light). Therefore, most of the 
electromagnetic energy flows out of the sample, leading 
to strong radiation losses. To overcome this disadvan- 
tage, new systems, e.g., layered superconductors instead 
of metals, should be considered in order to observe sur- 
face waves. 

Layered superconductors are strongly anisotropic me- 
dia where the current along the crystallographic ab- 
plane is similar to the current in bulk superconduc- 
tors, whereas transport along the c-axis is caused by 
the intrinsic Josephson effect through the layers (see, 
e.g., Refs. liliill ). Therefore, layered superconductors are 
characterized by two different plasma frequencies^ 7 -. One 
of them is of the order of oj p for normal metals, but the 
other one is much smaller and it coincides with the so- 
called Josephson plasma frequency uj ~ 1 THz. As was 
shown in Rcfs. [l8lri9l . the lower plasma frequency, ojj, 
defines the spectrum of surface waves in layered super- 
conductors. Namely, due to this feature the dispersion 
curves of the surface waves in layered superconductors, 
in contrast to normal metals, deviate substantially from 
the light-line, and these waves are well localized in the 
THz frequency range. Surface Josephson plasma waves 
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Figure 1: (Color online) Schematic geometry for the excita- 
tion of oblique surface waves. Here k*, k r , and k are the wave 
vectors of the incident, reflected, and oblique surface waves. 
Note that here S=superconductor and I=insulator. 



(SJPWs) can propagate along the interface between the 
external dielectric and layered superconductor for both 
geometries, when the superconducting layers are paral- 
lel to the sample surface or are perpendicular to it. As 
shown in Refs. [l9l - l2ll the spectrum of SJPWs propagat- 
ing along the layers consists of two branches, one above 
uij and the other below it. The spectrum of surface elec- 
tromagnetic waves propagating across the layers was pre- 
dicted and studied analytically in Refs. [2M23I 

In this paper, we study theoretically the excitation of 
the oblique surface waves (OSWs) propagating at some 
angle i9 with respect to the crystallographic c-axis, which 
is assumed to be parallel to the interface between a di- 
electric and a layered superconductor. We analyze the 
so-called attenuated-total-reflection (ATR) method for 
the OSW excitation using an additional dielectric prism 
with permittivity e p . The external electromagnetic wave 
from a dielectric prism is incident on a superconductor 
separated from the prism by a thin dielectric gap with 
permittivity Ed < e p (see Fig. [TJ) . In the absence of a su- 
perconductor, the incident wave completely reflects from 
the bottom of the prism, if the incident angle (p exceeds 
the limit angle tpt for total internal reflection. However, 
the evanescent wave penetrates under the prism a dis- 
tance about a wavelength, reaching the layered super- 
conductor and penetrating it. The wave vector of the 
evanescent mode is along the bottom surface of the prism 
and its value is higher than uj^Je2/c. This feature is the 
same as for surface waves. Thus, for a certain resonance 
incident angle, the spatial-and-temporal matching (the 
frequencies and wave vectors coincide) of the evanescent 
wave and the surface wave takes place. 

It is important to note that the excitation of the 
oblique surface wave is accompanied by another impor- 
tant resonance phenomenon. Evidently, the anisotropy of 
the current-carrying capability along the surface of the 
superconductor results in a conversion of the polarization 



of the terahertz radiation after reflection from the bound- 
ary of the layered superconductor. For example, if the 
incident wave has a transverse magnetic (TM) polariza- 
tion (a wave with the magnetic field parallel to the sam- 
ple surface), the reflected electromagnetic field contains 
a transverse electric (TE) wave (with the electric field 
parallel to the sample surface). Moreover, as was shown 
in Ref. l24l the complete transformation of the wave po- 
larization can be observed at an appropriate choice of the 
direction of the incident wave vector. 

We stress that this phenomenon is different from the 
Brewster effect. In the latter, light with TM polarization, 
for a definite incidence angle, is completely transmitted 
through a transparent dielectric surface, and, thus, the 
reflected light has only TE polarization. However, in the 
phenomenon presented here, the reflected wave has a TE 
polarization when the incoming bwave does not contain 
the wave with this polarization, i.e., when the incoming 
light is completely TM polarized. In other words, the 
reflected wave is TE-polarizcd because of the conversion 
of the incident TM wave, but not because of the separa- 
tion of the TE wave from the incoming mixed (TM+TE) 
wave. 

The ratio of the amplitudes of the reflected TE and 
TM waves is controlled by the wave frequency w, an- 
gles 1? and <p, and the thickness h of the dielectric gap. 
The main point is the to choose the angles 1? and <p for 
the incident wave such that the amplitude of the TM re- 
flected wave vanishes. It is known that, in the case of 
an isotropic metal, one can select the optimal value of 
the gap thickness h to provide the complete suppression 
of the reflected wave in resonance. Similarly, in the case 
of an anisotropic superconductor, we can select the best 
value of the angle $ to provide the complete suppression 
of the reflected TM wave in resonance. In the latter case, 
one part of the incident energy flux dissipates in the lay- 
ered superconductor and the other part is reflected as 
the TE wave. The theoretical study of the interplay of 
the OSWs excitation and the transformation of the wave 
polarization after its reflection from the boundary of the 
layered superconductor is the main goal of this paper. 

Thus, a layered superconductor plays two important 
roles in the phenomenon studied in this paper, which 
a simple metal or an isotropic superconductor cannot. 
First, due to the smallncss of the Josephson current 
across the layers, the excited oblique surface waves are 
well localized in the terahertz frequency range. Second, 
due to the strong current anisotropy in the surface of the 
layered superconductor, the excitation of oblique surface 
waves is accompanied by the transformation of the in- 
cident wave polarization after its reflection. Moreover, 
the complete transformation of the polarization can be 
achieved by the appropriate choice of the direction of the 
wave incidence. 

In the next sections, we derive and analyze the dis- 
persion relation for the oblique surface waves, and then 
study analytically the OSWs excitation by the ATR 
method and the transformation of the electromagnetic 
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Figure 2: (Color online) Geometry for the oblique surface 
waves. 



It is suitable to present the electromagnetic field as a 
sum of two terms that correspond to the ordinary and 
extraordinary evanescent waves, 



E(y)=^°>(y)+^ e \y), H(y) = H^> (y) +U^> (y) . (2) 

The electric field of the ordinary wave and the mag- 
netic field of the extraordinary wave arc orthogonal to 
the anisotropy z-axis: 



E<?\y) = Hf>{y) = Q 



(3) 



wave polarization. For this purpose, we calculate the 
electromagnetic field in the dielectric prism, in the dielec- 
tric gap, and in the layered superconductor, by solving 
the Maxwell equations with appropriate material equa- 
tions and boundary conditions. The electromagnetic field 
in the dielectric prism is presented as a sum of the inci- 
dent TM-polarized wave and two reflected waves with 
TM and TE polarizations. In the dielectric gap, we con- 
sider four evanescent waves, namely, increasing and de- 
creasing TM- and TE-polarized waves. The evanescent 
field in the layered superconductor consists of the ordi- 
nary and extraordinary decreasing waves. Then, using 
the continuity conditions for the tangential components 
of the electric and magnetic fields at the two interfaces 
(i.e., the dielectric gap-dielectric prism and the layered 
superconductor-dielectric gap), we express all wave am- 
plitudes via the amplitude of the incident wave. Using 
these expressions, we derive the reflection coefficient Rtm 
for the TM wave and the conversion coefficient Ttm-s-te 
for the mode transformation from TM to TE, as well as 
the expression for the portion A of the energy flux which 
comes to the excitation of the oblique surface wave and 
then dissipates in the superconductor. 



II. DISPERSION RELATION FOR THE 
OBLIQUE SURFACE WAVES 

Let the interface dividing an isotropic dielectric and 
a layered superconductor be located in the y = plane 
(see Fig. [2]) . We assume that the dielectric region (y < 0) 
and the layered superconductor region (y > 0) are semi- 
infinite along the y-axis. Both media are nonmagnetic 
ones. 

We choose the coordinate system in such a way that 
the z-axis coincides with the crystallographic c-axis of 
the superconductor. We study oblique surface waves, 
propagating along the interface dielectric-layered super- 
conductor at an angle i9 with respect to the c-axis, with 
the electric E and magnetic H fields in the form, 

E(x, y, z, t), H(x, y, z, t) oc cxp[i(np ~ ut)], (1) 

where p = (x, z) is the radius vector in the xz-plane and 
K = {k x ,k z ) is the wave vector of the oblique surface 
wave, |k| > k = lj/c. 



For each of these waves in the dielectric, the Maxwell 
equations give the relations between the components of 
the electric E^ and magnetic fields: 
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k~ dx ' 



(5) 



Here k = lj/c. In addition, the Maxwell equations give 
the same exponential law for the decay of both waves 
along the ?/-axis, 

E<£\y), E< e) (y), H { d °\y), H< e) (y) cx exp( l fc^), (6) 

with imaginary normal component of the wave vector, 

k dy = -i{k 2 x + k 2 z - Edk 2 ) 1 ' 2 . (7) 

The electromagnetic field inside the layered supercon- 
ductor is determined by the distribution of the gauge- 
invariant phase difference \ of the order parameter be- 
tween neighboring layers. This phase difference can be 
described by a set of coupled sine- Gordon equations (see, 
e.g., Refs. Il8ll25l - l30l) . In the continuum and linear ap- 
proximation, x can be excluded from the set of equations 
for the electromagnetic fields, and the electrodynamics 
of layered superconductors can be described in terms of 
an anisotropic frequency-dependent dielectric permittiv- 
ity with components e c and e a b across and along the lay- 
ers, respectively. Indeed, as was shown in Ref.LLZl the di- 
rect solution of the set of coupled sine-Gordon equations 
in the continuum and linear approximations is equivalent 
to the solution of the Maxwell equations for the uniform 
anisotropic medium with 



e c (fi) 



i 1 , ■ Vc 
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2 ^ab 

f2 



(8) 



Here we introduce the dimcnsionless parameters, fi = 
w/wj, ^ ab = Airaab/esUJjj 2 , and i/ c = 4*jr<x c /e a u}j; 
7 = A c /A 6 ^> 1 is the current-anisotropy parameter, 

1 /2 

A c = c/ojjEs and A a fc are the magnetic-field penetra- 
tion depths along and across the layers, respectively. 
The relaxation frequencies v a b and v c are proportional 
to the averaged quasi-particle conductivities a a b (along 
the layers) and a c (across the layers), respectively; ujj = 
(STreDjc/he) 1 / 2 is the Josephson plasma frequency. The 
latter is determined by the maximum Josephson current 
density j c , the interlayer dielectric constant s s , and the 
spatial period D of the layered structure. 

Contrary to the dielectric, the Maxwell equations give 
different laws for the decays of the ordinary and extraor- 
dinary waves in the layered superconductor, 

E(°)(y), H(°)( 2/ )(xexp(^(°)y), 



with 



E< e >(y), H< e >(i,)oc«p(i*#») (9) 

4? = *(*£ + fc " - k 2 s ab ) 1/2 , (10) 
k$ = i{k 2 x + k 2 z e c /e ab - k 2 e c fl 2 . (11) 



The relationships between the field components in the 
layered superconductor are as follows: 
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Figure 3: (Color online) Dispersion curves Q(ck/ljj) of the 
OSWs for e d = 4.56, e s = 16, 7 = 200, and v ab = u c = 
0. Curves 1-6 correspond to <d = 0°, 10°, 20°, 30°, 60°, 
and 90°, respectively. The hollow circles at Q. — fi cn d ~ 
fli = ^e s /{es — £d) show the end-points of the dispersion 
curves. In these points, the extraordinary wave, which is 
part of the OSW, is transformed into a bulk wave. When 
K — > 0, curve 6 as ymptotically approaches the frequency 
H = i^/e a /{e s +e d ) > 1. 



show the end-points of the corresponding dispersion 
curves. In these points, the extraordinary wave, which 
is present in the oblique surface mode, is transformed 
into a bulk wave. Indeed, the wave vector k^J vanishes at 
the end-points and is real and positive for 51 > fi on d. The 
analysis of Eq. (Q~4|) shows that fi C nd ~ ^1 for angles -d not 
very close to zero or ir/2. For very small d -C I/7 -C 1, 
we have O e nd ~ 1, and O e nd -> "f\/£ s /(e s + e d ) » 1 for 
(tt/2 - 0) « I/7. 

In the main approximation with respect to the param- 
eter 7> 1 and neglecting the dissipation, Eq. (fT4| gives 



Dn = 



1 



4y I 



0. 



(15) 



III. 



k x k z 



EXCITATION OF THE OBLIQUE SURFACE 
WAVES USING THE ATR METHOD 



rr(e) _ fcfc 4 We) rr(e) _ A e , r(<=) ^131 
J1 sx ~ 11 c ab J ^ S x i sy — i ^ab^sx \ LO I 

K x K z K z 

for the extraordinary waves. 

From the continuity conditions for the tangential com- 
ponents of the electric and magnetic fields at the interface 
y = 0, we derive the dispersion equation for the OSWs, 

klkil\£ ab -e d ) + e ab (k dy -k^)(k dy k^-kl) = 0. (14) 

Figure [3] illustrates the numerically calculated disper- 
sion curves of the OSWs for different values of the angle 
■d. The hollow circles at f2 = fiend ~ Oi = \J £s/(e s — e d ) 



In this section, we study analytically the OSWs exci- 
tation using the ATR method. Let the structure under 
study be comprised of a dielectric prism with permittivity 
e p , a dielectric gap with thickness h and permittivity e d , 
and a layered superconductor described by a frequency- 
dependent diagonal tensor of dielectric permittivity with 
components e xx = s yy = e ab and e zz = e c given by 
Eq. ([EJ (see Fig. [I]). All media are supposed to be non- 
magnetic. We define the coordinate system so that the 
prism region occupies the half-space y < 0, the super- 
conductor region occupies the half-space y > h, and the 
z-axis coincides with the crystallographic c-axis of the 
layered superconductor. 
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Consider now a TM polarized wave incident from a di- 
electric prism on a layered superconductor. The incident 
angle ip exceeds the limit angle ipt = arcsin[(£ ( j/£ p ) 1 / 2 )] 
for the total internal reflection. In this case the incident 
wave can excite a surface wave if the resonance condi- 
tion. This condition implies that the tangential compo- 

1/2 

nent ke p shop of the incident wave vector is equal to 
the wave vector k of the oblique surface wave, 



he 1 ' 2 



Here and henceforward, we omit the multiplier 
exp[i(np — uit)] in all expressions for the electromagnetic 
field components. 

For the dielectric gap region, we define the electro- 
magnetic field in the same manner as for the prism re- 
gion, namely, as the superposition of the TM- and TE- 
polarizcd waves with the following components: 



E 



dx TM 



Bi exp(ik dy y) + B 2 exp(-ik dy y), (22) 



The electromagnetic field in the dielectric prism is a 
sum of three terms that correspond to the incident TM 
polarized wave and two reflected waves with the TM and 
TE polarizations. Thus, the cc-component of the electric 
field in the prism can be presented as 

E px (x,y,z,t) = [E x exp(ik py y) 

+ (Ktm + E r xTE )exp(-ik py y)]exp[i{Kp- ut)] (16) 
where E x , E r x 

tm ! an d E r x TE are the amplitudes of the 
x-components of the electric field for the incident waves 
and for the reflected TM and TE waves, respectively; 
kpy = ke p ^ 2 cos ip is the normal component of the wave 
vector of the incident wave. 

The electric and magnetic field components in the 
prism can be expressed via E x , E xTM , and E r x TE using 
the Maxwell equations, 



E — 
j^ py — 



„2 



hpy k x 



E l x cxp(ik py y) - E r x TM exp(-ik py y) 



(17) 
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P z - T~E X exp(ik py y) 

K x 



~~r~E T 

yk x 



TM 



-r E l TE ) exp(-ik py y), (18) 



= p . z E x exp(ik py y) 



kx k' 
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£ p kk z 
~kl~ x 



TM 



h h 

kk~ 



Hpy = -^r Er xTE cxp(-ik py y), 



e Ite) exp(-ifc py y), (19) 



(20) 



H,, 



'l~ E x ex P( ik pyy) 

K py 



' ke 



-El 



TM 



-El 



TE 



y i>y 



cxp(-ik py y). 



(21) 



E 



d y TM 



kd y k x 



Bi exp(ik dy y) - B 2 exp(-ik dy y) 



Edz TM — ~j~~E dx TM, 



(23) 
(24) 



e d kk z 



dx TM 



— E dy TM, 11 dz TM = J — £j d y TM, 



H, 



s d kk x 



(25) 



E dx te = C\ exp(ik dy y) + C 2 exp(-ik dy y), (26) 



E dz 



TE 



-E dx TE, 



(27) 



H, 



dx TE 



kx k dy 
kk z 



Ci cxp(ik dy y) - C 2 exp(~ik dy y) 



(28) 



K k z 

H dy TE = T-r~H dx TE, H dz TE — ~~r~H dx TE- (29) 

The electromagnetic field in the layered superconduc- 
tor is the superposition of ordinary and extraordinary 
waves, and it is described by the same formulaes as in 
the previous section [see Eqs. ([<?]) - (flU)) ]. 

Using the continuity boundary conditions for the tan- 
gential components of the electric and magnetic fields at 
the two interfaces (i.e., the dielectric prism-dielectric gap 
and the dielectric gap-layered superconductor) , we obtain 
eight linear algebraic equations for eight unknown ampli- 
tudes (for 4 waves in the dielectric gap, 2 waves in the 
layered superconductor, and amplitudes £J TM , -E^te of 
the reflected waves in the dielectric prism) . Solving these 
equations, we derive the reflectivity coefficient 



#tm (</=>)■#) 



|E r TM (y,fl)l s 



for the TM wave, the coefficient 

|E r TE (^)| 2 



Ttm->te(v,'') 



' TEV 

IE 



'■12 



(30) 



(31) 



of the transformation of the incident TM wave to the 
reflected TE wave, and the absorptivity 



A(p,r9) = l-R TM (p,tf)-T 



TM-vTE 



{#,■&). (32) 
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It is suitable to present these coefficients in the form, 
(Do + Dc) 2 + (Ltm - L TE - T) 2 



R 



TM 



(D + D c 



J TM_ 
TM 



L 



J TE_ 
TE 



r) 



2 ' 



Ttm- 



4LtmLtm 



■TE 



(Do + D c ) 2 + (L TM + L TE + T) 2 ' 



A = 



atl 



TM 



(D + D c ) 2 + (L TM + L TE + T) 2 ' 



(33) 



(34) 



(35) 



The expressions for the parameters D c , Ltm, Lte> and 
r are very cumbersome. However, they can be signif- 
icantly simplified for the case of small dissipation and 
weak coupling of the electromagnetic fields in the lay- 
ered superconductor and the dielectric prism, when the 
following inequalities are satisfied: 



r< 1, 



C h = exp(-2|fc d j / |/i) < 1. 



(36) 



For this case, using the smallncss of the parameter 7 
we obtained the following formulaes: 



Dr. 



2C 4 cos 2 ■& 
£ P - £d 



\k d y\ 2 k 2 py + E ? d k 2 \k dy \ 2 + ^f(l 



I fcdy I 



^ + \k dy \ 2 k 2 



mi 



L 



TM 



^TE 



e d cos 2 e p e d k 2 1 k dy I k py 



£p £d K ^ 4" \kdy\ 2 k 2 y 



, (37) 



(38) 



2 a 



£d cos^ 1? k% k I I kpy 

Sp Ed k^ ~h \kdy\^kpy 



abt 



(39) 



Rtm(<P,$) occurs when the tangential component k of 
the incident wave vector coincides with the wave vector 
of the oblique surface wave. Indeed, the equation 



D n 



D r = 



(41) 



is the dispersion relation for the OSWs in the system 
under consideration, which consists of the layered super- 
conductor, the dielectric gap, and the dielectric prism. 
Because of the coupling of the electromagnetic fields in 
the superconductor and the dielectric prism, this relation 
differs slightly from the dispersion equation Do = for 
the OSWs in the system without the prism. 

The parameters Ltm, Ltej and T in the denominator 
in Eq. (|30|) define the resonance width. These are related 
to three channels for the losses of the electromagnetic en- 
ergy of the OSWs. The term T describes the dissipation, 
the term Ltm is related to the energy leakage through 
the dielectric prism in the form of a TM bulk wave, and 
the term Lte describes the leakage in the form of a TE 
wave. 

The resonance suppression of Ltm (<£>,$) can be com- 
plete if both terms in the numerator in Eq. (|33[) vanish 
at the same angles ip and -d. Namely, the complete sup- 
pression of Ltm^j^) occurs if the dispersion relation 
Eq. (|4"Tj) and the condition 



Ltm — L 



TE 



r. 



(42) 



are fulfilled simultaneously. Recall that, in the isotropic 
case, the condition for complete resonance suppression of 
Rtm can be realized for metals or superconductors 

with finite dissipation parameter T (see, e.g., Ref. 0). 
However, in the anisotropic case considered here, this 
phenomenon can be observed even when T = 0. In the 
main approximation with respect to the small parameters 
7 _1 <C 1 and Ch <C 1, the corresponding angles ip and i9 
can be found by solving the set of equations which follow 
from the equalities Do = and Ltm = Lte, 

• 2 2 -2 £ d 
sm if cos ip = sin ip 

Ep 



sm ip 



< sm ip sm ip 



E,n 2 



o 2 



(43) 



e* \k dy \ r _ n 2 k 2 x {\k dy \ + \k% ] \) 
e 2fi|fcW|' Lab ~ Vab 2 1 k 3 

(40) 



IV. RESONANCE SUPPRESSION OF THE 
REFLECTIVITY 

Equation ([33]) together with Eqs. (I5J), (|57 )l -(|3D )l de- 
scribes the resonance behavior of the reflectivity coeffi- 
cient -Rtm (</?,$) when changing the direction of the in- 
cident wave propagation. The resonance suppression of 



o 4.2 a 4" Ed . 4 2 Ep + Ed . 2 -, / a a\ 

2 cot § = — — sm ip -\ sin <p — l. (44) 

£d £d £d 

The graphical solutions of these equations are presented 
in Fig. U for different material parameters. Figure [5] 
demonstrates the complete resonance suppression of the 
reflection coefficient that occurs for the dissipationlcss 
case, namely at the angles tp = <p 4 m t and d — i?mt, which 
are the solutions of the set of equations (j4"3"|) , ([4"4"]l shown 
in Fig. |U 

For the lossy superconductor with r ^ 0, the complete 
resonance suppression of the reflectivity can also be ob- 
served. In this case, the energy flux of the incident TM 



7 




50 
0>(deg) 



Figure 4: (Color online) Graphical solutions of the set of equa- 
tions (|43[) and (|44[) . The red solid curves correspond to the 
dispersion relation for OSWs and the blue dashed curves are 
for the condition Ltm = Lte- The curves start from the 
angles tp t , which are the critical angles for the total inter- 
nal reflection from the bottom of the dielectric prism. The 
material parameters are chosen as e p = 11.6 (this value cor- 
responds to silicon, which is frequently used as the material 
for the dielectric prism for the ATR experiments in the THz 
frequency ranged), £ B = 16 (this interlayer permittivity corre- 
sponds to Bi2Sr2CaCu208+,r^). Curves are plotted for Ed = 1 
(the vacuum gap, curves 1), Ed = 2.04 (teflon, curves 2), and 
Ed = 4.56 (quartz, curves 3). Other parameters are Q. = 1, 
v a b = = 0. The intersections of red and blue curves occur 
at tp = <^>j Ilt ~ 22.5°, $ = $i n t « 30.6° (curves labeled by 1); 
pint « 34.9°, tfint w 28° (curves 2); tp int « 76.4°, # int « 21.6° 
(curves 3). 
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wave is partially reflected as the TE wave. The other part 
A of the energy [absorptivity A in Eq. (|32|)] is dissipated 
in the superconductor (see Fig. [B]). 

It is important to note that the phenomenon described 
here owe the same origin as the wave interaction with an 
open resonator—. The tunneling of an incident plane 
wave through an open ID resonator is characterized 
by the reflection and transmission coefficients R and T 
which play the same role as the coefficients Rtm( ( P>'&) 
and Ttm->te in our analysis. Moreover, the formulas for 
the coefficients R and T presented in Ref. [3l| are exactly 
the same as Eqs. ([33]) and ([34]). The role of the leakage 
parameters Ltm and Lte are played by the inverse leak- 
age Q-factors which are related to the transmittances of 
the two tunnel barriers discussed in Ref. [3lT However, 
it is necessary to emphasize a principal novel feature of 
the phenomenon considered here. We now pay attention 
to the dual role of the coupling of the electromagnetic 
fields in the dielectric prism and in the layered supercon- 
ductor in the conversion of the terahertz wave polariza- 
tion. First, due to the coupling, the incident TM wave 
resonantly excites the oblique surface wave. Then, the 
OSW supplies the TE bulk wave by the energy which 



Figure 5: (Color online) Complete resonance suppression of 
the reflectivity in the dissipationless regime. The depen- 
dence of the reflectivity 7?tm and the transformation coef- 
ficient Ttm-j-te on the incident angle ip are shown by the 
solid red and dashed blue curves for the same parameters as 
in Fig.g] (a) e d = 1, ■& = 30.6°; (b) £ d = 2.04, ■& = 28°; (c) 
Ed = 4.56, •& = 21.6°. These three cases correspond to the 
three intersections in Fig. f3] The gap thickness h is equal to 
1.2c/wj. 

passes through the dielectric gap due to the same cou- 
pling. Thus, the oblique surface wave serves as an inter- 
mediary between the incident TM and the reflected TE 
waves in the phenomenon of the transformation of the 
wave polarization. 

V. CONCLUSION 

We have theoretically studied oblique surface waves 
which propagate along the interface between a dielectric 
and a layered superconductor. We consider a geometry 
in which the interface is perpendicular to the supercon- 
ducting layers. We have analyzed the case of arbitrary 
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Figure 6: (Color online) Dependence of the reflectivity -Rtm 
(red solid curves), the transformation coefficient Ttm->tb 
(green solid curves), and the absorptivity A (dashed blue 
curves) on the incident angle ip for lossy superconductors for 
Q = 1, e s = 16, 7 = 200, h = 0.5c/uj, e p = 11.6, e d = 1. The 
values of other parameters are: (a) v a b = v c = 10 -5 , i? = 6.3°; 



(b) Vab 



10" 



i? = 26° 



direction of the OSW propagation with respect to the 
layers. In this general case, the electromagnetic field of 



the OSWs is a superposition of ordinary and extraordi- 
nary waves. We derived the dispersion equation for the 
OSWs, which is the most general dispersion relation for 
surface Josephson plasma waves in the geometry consid- 
ered. We have also shown that the dispersion curves 
have end-points where the extraordinary mode trans- 
forms from evanescent wave to bulk wave propagating 
deep into the superconductor. In addition, we have an- 
alytically solved the problem of the resonance excitation 
of the OSWs by the attenuated-total-reflection method 
using an additional dielectric prism. Due to the strong 
current anisotropy in the boundary of the superconduc- 
tor, the excitation of the OSWs is accompanied by an 
additional important phenomenon. Namely, the electro- 
magnetic field with orthogonal polarization appears in 
the wave reflected from the bottom of the prism. We 
show that, for definite optimal combinations of the prob- 
lem parameters (the wave frequency, the direction of the 
incident wave vector, the thickness of the gap between 
the dielectric prism and superconductor, etc.), there is a 
complete suppression of the reflected wave with the same 
polarization as the incident wave. This phenomenon is 
an analog of the so-called Wood anomalies in the re- 
flectance of isotropic metals known in optics. However, 
in the isotropic case, the optimal thickness of the dielec- 
tric gap between the prism and the metal sample should 
be chosen for the complete suppression of the reflected 
wave. However, in the anisotropic case considered here, 
this phenomenon can be more easily observed by the ap- 
propriate choice of the direction of the wave propagation 
with respect to the superconducting layers. We obtained 
the conditions for the observation of the complete sup- 
pression of the reflectivity and found that, contrary to 
the isotropic case, this phenomenon can be observed even 
in the dissipationless limit. In such a regime, a complete 
conversion of the wave polarization takes place, i.e., 100% 
of the energy flux of the incident TM wave reflects from 
the bottom of the prism as a wave with TE polarization. 
Similar conversion phenomena can be observed for the in- 
verse transition, from TE to TM modes, as well as for the 
transition from incident ordinary waves to extraordinary 
ones (or viceversa). 
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